Abstract 

We analyze electrostatic interaction between a 
sharp conducting tip and a thin one-dimensional 
wire, e.g., a carbon nanotube, in a scanned gate 
microscopy (SGM) experiment. The problem is 
analytically tractable if the wire resides on a thin 
dielectric substrate above a metallic backgate. 
The characteristic spatial scale of the electro- 
static coupling to the tip is equal to its height 
above the substrate. Numerical simulations in- 
dicate that imaging of individual electrons by 
SGM is possible once the mean electron sepa- 
ration exceeds this scale (typically, a few tens of 
nm). Differences between weakly and strongly 
invasive SGM regimes are pointed out. 
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Scanned gate microscopy (SGM) is a recent 
addition to the arsenal of modern scanned probe 
techniques that enable one to study and ma- 
nipulate nanoscale objects at a single-electron 
level. Much interest has been attracted by exper- 
iments pmisiiiiisjin] where the SGM has been 

employed for probing carbon nanotube (CNT) 
quantum dots under the conditions of Coulomb 
blockade. In those experiments, a movable sharp 
tip with a controllable electrostatic potential was 
employed to add or remove electrons from the 
dots one by one, which was monitored via trans- 
port measurements. 

The goal of this work is to explore the capa- 
bilities of the SGM for imaging the real-space 
electron structure in a CNT dot. In parallel, 
we develop an analytical approach that clarifies 
what controls the fundamental limits of the SGM 
spatial resolution. We focus on a typical SGM 
geometry where the CNT resides on a dielectric 
substrate above a metallic backgate (Fig. . We 
show below that if the substrate is thin, so that 
the gate efficiently screens Coulomb interactions, 
the requisite electrostatic problem is tractable. 
We find the potential created by the charged tip 
in terms of elementary functions. For realistic 
experimental parameters our formula proves to 
be accurate to about 5%, which enables us to test 



with confidence the assumptions made about the 
same quantity in prior literature [U Ej . In par- 
ticular, we show that for the case of a sharp tip 
the spatial resolution of the SGM is set by the 
tip's height above the substrate. For thick sub- 
strates, the resolution degrades and becomes of 
the order of the geometric mean of the tip height 
and the CNT length. Typical experimental pa- 
rameters put one in between these two limiting 
cases, with the net resolution of the order of a 
few tens of nm. We propose that this may be 
sufficient for imaging of individual charges in a 
few-electron quantum dots and support this idea 
by numerical simulations. Our method of anal- 
ysis can be extended to other one-dimensional 
(ID) and quasi-lD systems (multiple quantum 
dots, CNT networks, etc.), and so we hope it 
may be a useful tool in the experimental prac- 
tice of SGM. 

We consider a CNT of radius R resting on an 
insulating substrate of thickness D deposited on 
a metallic gate (Fig. We assume that the 
contact resistances are high, so that the CNT be- 
haves as a quantum dot. Such a dot is probed by 
a movable conducting tip, hovering at the height 
h above the point r t = (xt,yt) of the substrate 
plane. The radius w(() of the tip function 
of distance Q = z — h from its end is given by 
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Figure 1: (Color online) The SGM setup studied 
in this paper: a movable metallic tip that probes 
a ID wire (the horizontal cylinder). 

w(C) = 8( at I < C < L, where 8 < 1 is the 
apex angle and L is the full length of the tip. 
Parameter I <C /i accounts for the rounding of 
the tip's end 

Our first task is to analyze Coulomb inter- 
actions among the charges on the CNT, gate, 
and the tip. For simplicity, we set all dielec- 
tric constants to unity j^j. We start with a 
qualitative picture. The excess charge density 
Pd(x) induced on the CNT by the external gates 
can be deduced from classical electrostatics of a 
thin metallic wire. Earlier microscopic calcula- 
tions [7| ^2 El proved that this is a good ap- 
proximation when lengthscales of interest exceed 
e/pd. If the distance D to the screening gate is 
much smaller than the length of the CNT, pd{x) 
is nearly uniform. This is because the effective 
interaction potential 



U(x) 



\x\+R Vx 2 + AD 2 



(1) 



between electrons on the dot is short-range. For 
R <C D the capacitance per unit length of the 



CNT is c d = 1/2 In (2D/R). Actually, even with- 
out the screening gate, i.e., at D = oo, pd(x) 
is nearly uniform, except near the ends |13| . 
because the potential U(x) oc is only 

marginally long-range in ID. 

Let us now discuss the linear charge density 
Pt(z) of the tip. Interestingly, it also becomes 
uniform ^1] if h and D are both small and the 
tip is sharp, 8 <C 1. In this case w(Q is a 
slowly varying function; therefore, a segment of 
the tip at height z above the substrate is sim- 
ilar to a piece of thin wire of radius w(z — h). 
The distance to the screening gate for this seg- 
ment is z + D. A quick estimate of the ca- 
pacitance of the tip per unit length q can be 
obtained from the formula for c d by replacing 
R — > w(z — h) ~ [z — h)8 and D — ► z + D, i.e., 



2 ln[2(z + D)/w(z - h)] 2 ln(2/0) 



• (2) 



This result is confirmed by a formal perturbation 
theory that yields pt(z) as a power series of the 
small parameter 1/| ln#|, similar to Ref. ^31- in 
the leading order we obtain (£ = z — h) 



1 

Pt 



In 



4(£ ~ QC C + 2h + 2£> 
L + ( + 2h + 2D w 2 (() 



■ (3) 



This formula applies at all £ except very near the 
ends of the tip, £ = 0, L. We tested it by solv- 
ing the electrostatic problem for the tip and the 
gate numerically. An excellent agreement was 
reached for, e.g., 8 = 0.1 ~ 6° and h = 30 nm, 
which are not too difficult to achieve experimen- 
tally (HI. For such 8 and h the nonuniformity 
of pt within the important region < £ < 2D 
(see below) is about 15% for D = 30-200 nm. 
These are the parameters adoped in our subse- 
quent simulations described shortly below. 

The approximate uniformity of pt entails a 
simple model for the potential Ut(r) induced by 
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the tip on the substrate. This potential is cre- 
ated both by the tip itself and by its oppositely 
charged image [Fig. 12(a)] • It is convenient to 
move the latter into the z > half-space by re- 
flection [Fig. |2[b)]. As a result, the two charge 
densities cancel each other almost everywhere. 
The important uncompensated piece is a uni- 
formly charged rod of length 2D [Fig. [2(c)]. The 
potential induced by this rod is 



U t = ctVtl sinh 



h + 2D 
Ar 



h 

Ar 



sinh- 1 — , (4) 



where Ar = | r — | is the in-plane distance from 
the tip and Vt is the tips's voltage. In partic- 
ular, for a thin dielectric substrate, Ut has the 
Coulomb form, 



U t (r) = Q/^A^ + h 2 , D^h. (5) 

In this limit the effect of tip reduces to that of 
a point charge Q = DVt/ln(2/9) positioned at 
height h above the substrate [Fig. 12(c)] . Clearly, 
h sets the range of the tip-dot interaction and is 
thus the sole geometric parameter that controls 
the SGM spatial resolution (Res). 

As D increases, the resolution degrades. The 
following estimate can be derived: 



(6) 



Res ~ mm{h + D, Vbh}, 



where b 3> h is the half-length of the quantum 
dot probed by the SGM. The reason why b enters 
at large D is the slow decay of Ut with Ar. In 
this case Res is determined not by the behavior 
of Ut itself but by that of its variation across the 
system, AC7 t (r) = Ut{r) — Ut(b,0). Indeed, an 
additive constant in the dot-tip interaction en- 
ergy has no effect on spatial resolution. A rough 
estimate of Res is then the value of Ar at which 
AUt drops to a half of its maximum value at 
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Figure 2: Geometrical construction that enables 
one to replace a long tip by a uniformly charged 
rod of a much smaller length 2D (see main text). 



Ar = 0. Thus, in the case of a thick substrate, 
D 2> b, where Eq. (jlj) yields 



AU t ^ c t V t In 



VAr 2 + h 2 + h 



r t < b, 



(7) 

the second expression in Eq. © is recovered. 

Note that the actual resolution of the SGM 
depends on temperature. At low T Coulomb 
blockade spectroscopy enables one to detect very 
small variations of the energy of the system as 
a function of r^. Hence, the spatial resolution 
also improves. Nevertheless, outperforming the 
estimate © by more than a factor of two or so 
would probably require impractically low T. 

The substrates used in the current SGM ex- 
periments commonly have the effective thickness 
of D ~ 50 nm [S] , which is not much larger than 
a typical h. Therefore, Eq. © is acceptable as a 
first approximation. This also explains why an 
analogous model for the tip — a small metallic 
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sphere — postulated previously was in a 
reasonable agreement with the experiment. 

Armed with Eq. @ , we now turn to the com- 
putation of representative SGM patterns. The 
conductance measured in an SGM experiment is 
directly related to the tip-dot capacitance Cdt- 
When the preferred dot charge Qd = CdtVt + 
const is a half- integer multiple of e, the Coulomb 
blockade is lifted and the conductance peak is 
measured j!5j . The problem reduces to calcula- 
tion of Cdt(rt). Let us first discuss the case of a 
metallic CNT where it can be done analytically. 

The key idea is to treat the inverse self- 
capacitances of the dot and the tip, \n{2D/R) 
and ln(2/6>), respectively, as large parameters. 
Then Cdt can be computed by a perturbation 
theory. In the leading order one needs to inte- 
grate the product Ut(x)pd(x) over the length of 
the dot. This gives the tip-dot interaction en- 
ergy from which Cdt can be extracted. In the 
simplest case D <C h one can use Eq. © to get 

DlQ f xt+b+yj '{xt+bf +yZ+h? \ 
r \x t -b+y/(x t -b?+yl+b? ) 

dt ~ 2\n(2/9)\u{2D/R) lj 

In deriving this equation, additional screening 
by source and drain leads was neglected. This 
is legitimate away from the leads, at b — \xt\ 3> 
D,hi, where hi is the height by which the leads 
rise above the substrate. 

The contours of Q d (r t ) = C dt V t = (2V + 
l/2)e = const computed according to Eq. (jHJ) for 
typical experimental parameters are presented in 
Fig. |21 They are oval-shaped near the CNT and 
become nearly circular at larger N, as in exper- 
iment f3j. To verify our results for Cdt quantita- 
tively, we compared them with the output of the 
numerically exact capacitance calculator FAST- 
CAP [TS|. For parameters specified in Fig. and 
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Figure 3: High-conductance contours CdtVt = 
(N + 1/2) e = const as function of the tip posi- 
tion for a metallic CNT. N runs through a set of 
consequtive integers. Parameters: b = 250 nm, 
R = 2 nm, D = 30 nm, and V t = 0.3 V. 

the tip of length L = 10 pm we found a 5% agree- 
ment, which is gratifying given the simplicity of 
our approach and the approximations made. 

Classical electrostatics we used thus far is valid 
for a metallic nanotube where the average inter- 
electron separation a is small. A more interest- 
ing question is whether the SGM can help us to 
see effects beyond continuum electrostatics, for 
example: Can SGM image individual electrons? 

As explained earlier, Res is of the order of the 
height h of the tip above the substrate, typically, 
a few tens of nm. This is insufficient for resolv- 
ing individual electrons in a metallic CNT. On 
the other hand, in a semiconducting CNT the 
mean inter-particle separation a = 20-30 nm can 
be achieved. In such CNTs electrons (or holes) 
can be added to the conduction (valence) band 
one by one starting from zero ^7j. Since the 
typical size of a CNT quantum dot is a frac- 
tion of a micron, it would contain only a few 
charge carriers altogether. In this regime the 
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electrons form a Wigner molecule: a state where 
they are highly localized near the classical equi- 
librium positions. The Wigner-molecule approx- 
imation is accurate if r s = a/2as > 4, [T%1 IT§] 
where as is the effective Bohr radius. At large 
r s the charge distribution of each electron can be 
approximated by a Gaussian with the standard 
deviation w « 0.52aB(a/as) 3//4 . For example, 
for a ~ 40 nm (as in Fig. @J and as ~ 1.4 nm 
(as in the experiment ^7J) we get w ~ 6.4 ag, in 
agreement with numerical simulations of Wigner 
molecules [TBI. 

Suppose Res < a, then as the tip moves along 
the dot, it interacts primarily with the nearest 
electron, and so continuum electrostatics is no 
longer valid. Instead, one expects a significant 
oscillatory modulation of the tip-dot coupling: 
maxima when the (repulsive) tip is positioned 
directly above an electron and minima when it is 
between two adjacent electrons. Such variations 
may be detectable by examining the spacing of 
the conductance peaks. 

To verify these qualitative ideas we carried out 
a series of numerical simulations. The Hamilto- 
nian of the system H was taken to be 

H= ^2 U(xj - Xj) + ^ eUt(xj) + H COTi , 



N 



H C on = y^e[0 con (xt) + V g ), 



(9) 



i=i 



where electron-electron and electron-tip interac- 
tions are given by Eqs. Q and @, V g is the 
gate voltage relative to the that of the leads, 
and ^con(aj) is the confinement potential due to 
contact potential difference between the leads 
and the CNT. Following Ref. [IT], we chose 
0con(aO = {k/2)(x 2 — b 2 ), with the coefficient 
k close to the number estimated in that paper. 
The ground-state energy En and the electron 
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Figure 4: (Color online) Results of numerical 
simulations for SGM of Wigner molecules, (a) 
Conductance [estimated per Eq. ()1U[)] (b) chem- 
ical potential, and (c) ground-state electron po- 
sitions for h = 30 nm, D = 200 nm, and Vj = 
— 0.1V. In plot (c) every other curve from the 
top is for the six-electron molecule, the rest are 
for the five-electron state. At a higher tip volt- 
age, Vt = —0.25 V, plot (c) evolves into (f) 
and plot (b) into (e). At a larger tip height, 
h = 90 nm, graph (b) evolves into (d). Dia- 
monds correspond to x% at which the tip is di- 
rectly above one of the electrons. 
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positions Xj were found numerically for a dense 
grid of tip coordinates on the interval < xt < b. 
Figure 21 shows the results. Unlike Fig. Inhere we 
focus on line scans where the tip stays directly 
above the CNT, i.e., at y t = 0. 

At low T there are only two important charge 
states to consider, Qd = Ne and Qd = (N + l)e. 
In the simulations the gate voltage V g was ad- 
justed to have iV = 5. The dependence of the 
electron positions on Xt in the two competing 
configurations is shown in Fig. H(c) for the rela- 
tively low voltage of Vt = — 0.1 V. The displace- 
ments of the electrons by the tip are rather small, 
so here the tip is a weakly invasive probe. Shown 
by diamonds are the tip positions where it is di- 
rectly above an electron. The chemical potential 
\i = En+i — En for the same SGM scan is plot- 
ted in Fig. Efb). It exhibits the expected spatial 
variations: maxima and minima near x^'s of the 
diamonds. At the intersection of the n{xt) curve 
with the horizontal line /i = the two charged 
states become degenerate. At such points the 
conductance G through the dot has peaks. The 
calculation of the shape of these peaks is beyond 
the scope of this work. However, for illustrative 
purposes, in Fig. Ufa), we plot the expression 



2Ti? sd sinh(/x/T)' v ' 

which holds for large metallic dots |15| . Param- 
eter R s d here denotes the sum of the tunneling 
resistances at the source and the drain leads. 

According to the discussion above, imaging 
the electron positions amounts to finding the 
maxima and minima of the chemical potential. 
To achieve that in experiment one would need 
to repeat SGM scans at different V g . Let fJ-(xt) 
be the chemical potential at some initial gate 
voltage V*. At a different V g , the chemical po- 
tential changes by a constant, i.e., the curve in 



Fig. E^b) shifts up or down as a whole. The 
conductance peaks are found at the roots of the 
equation fi(x t ) = eV* — eV g . Based on this rela- 
tion, the entire curve of n{xt) can be deduced by 
tracking the peak positions as a function of V g . 

As the tip's voltage increases, it perturbs the 
system stronger, see Fig. ^f). Here we focus on 
the case of a repulsive tip, which expels electrons 
from underneath itself. This creates a double 
quantum dot: Dot 1 in front of the tip and Dot 
2 behind it [regions above and below the dashed 
line in Fig.0Jd), respectively]. As the tip moves, 
the following sudden changes in the double-dot 
system take place: (1) an electron is expelled 
from Dot 1 into the lead, which occurs when ji < 
changes to \x > 0, (2) another electron enters 
Dot 2 from the opposite lead, at the points where 
\i > changes to \i < 0, and (3) an electron 
jumps between Dots 1 and 2, at the positions 
labelled by diamonds in Fig. Eff). Under the 
assumptions that contact resistance R s d is higher 
than the inter-dot tunneling resistance and the 
temperature is low enough, the peaks in total 
conductance G coincide |2U] with events (1) and 
(2), i.e., the points fi = 0. 

The strongly invasive SGM can be easily dis- 
tinguished from a weakly invasive one experi- 
mentally. First, in the invasive regime n(xt) 
curve has cusps [Fig. HJe)] whereas in the weakly 
invasive one it varies smoothly [Fig. Efb)]. Sec- 
ond, in the double-dot case the peak heights of 
the conductance have a strong modulation (beat- 
ing pattern) |5J as a function of xt and V g . The 
highest G is achieved if the Coulomb blockade in 
the two dots is lifted simultaneously. This can 
be realized by tuning V g to some special values 
at which cusps of fJ,(xt) touch the \i = line. 
We do not attempt to discuss G in any more de- 
tail because its calculation for the double-dot is 
a complicated problem [201. Equation (|1U|) be- 
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comes a poor approximation in this regime. 

The effect of the tip height on spatial resolu- 
tion is demonstrated in Fig. Efd) : for h = 90 nm 
the oscillations in fi(xt) are completely obliter- 
ated. This confirms that the resolution limit 
strongly depends on h. In contrast, having a 
rather large D = 200 nm does not degrade the 
resolution much. It may be surprising at first 
but, in fact, Eq. © gives Res ~ 55 nm (us- 
ing b = max\xj\ ~ 100 nm), which is only 
slightly larger than the mean electron spacing 
a ~ 40 nm. Note also that a rather low temper- 
ature, T = 0.5 K, is chosen. 

Finally, let us discuss the results for a two- 
dimensional SGM scan, Fig. 0D We once again 
see contours of G, akin to Fig.|3]but with striking 
differences. One of the contours contains pro- 
nounced oscillatory features — "wiggles". An- 
other contour (closest to the CNT) has been bro- 
ken into three disconnected pieces, as though the 
amplitude of the corresponding wiggles exceeded 
this contour's original width. Similar to the ID 
scans of Fig. 21 such effects originate from dis- 
creteness of the electron charge that causes mod- 
ulation of the chemical potential as the tip passes 
by individual electrons. In general, the contour 
with N wiggles separates the interior region of 
N— 1 electrons from the exterior one with N elec- 
trons in the ground-state. The lack of wiggles on 
more distant contours is explained by noting that 
the spatial resolution is now set by \Jh 2 + y\ , 
i.e., it worsens as the tip is moved laterally away 
from the dot. At large yt continuum electrostat- 
ics applies, and so the distant G-contours are 
similar to those of Fig. G3 

To assess how stable the predicted pattern is 
against thermal smearing, G in Fig. is com- 
puted for T = 15 K. For simplicity, the calcu- 
lation is done based on Eq. (|10|) although this 
equation is a poor approximation in this case. 
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Figure 5: (Color online) Contours of high con- 
ductance in a low-T two-dimensional SGM scan. 
Simulation parameters: D = 30 nm, h = 30 nm, 
V t = -0.4 V, T = 15 K. 

In particular, it gives a constant intensity along 
the contours of high G whereas in reality it is ex- 
pected to have strong variations. A more accu- 
rate calculation of G is deferred for future work. 

In conclusion, we have studied electrostatics 
of a prototypical SGM experiment in which a 
CNT quantum dot is probed by a sharp metal- 
lic tip. We developed an analytical model for 
the coupling between the tip and charges on the 
sample plane and verified it by numerical simula- 
tions. Our model explains why the resolution is 
controlled primarily by the tip-dot distance (the 
height above the substrate and the lateral dis- 
tance away from the CNT) rather than by much 
larger overall dimensions of the tip and the dot. 
We predicted that the SGM of a semiconducting 
CNT in the regime of strong electron-electron 
interaction may be able to detect positions of 
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individual electrons forming a Wigner molecule. 
Lastly, the tip can act as either a weakly or a 
strongly invasive probe, and so the SGM can be 
employed both for single-electron imaging and 
for nano-manipulation. 

This work is supported by UCSD ASCR and 
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P. McEuen, and J. Zhu for valuable comments. 

References 

[1] S. J. Tans and C. Dekker, Nature (London) 
404, 834 (2000). 

[2] A. Bachtold, M. S. Fuhrer, S. Plyasunov, 
M. Forero, E. H. Anderson, A. Zettl, and 
P. L. McEuen, Phys. Rev. Lett. 84, 6082 
(2000). 

[3] M. Bockrath, W. Liang, D. Bozovic, 
J. H. Hafner, C. M. Lieber, M. Tinkham, 
and H. Park, Science 291, 283 (2000). 

[4] M. Woodside, PhD thesis, University of 
California, Berkeley, 2001; M. T. Woodside 
and P. L. McEuen, Nature (London) 296, 
1098 (2002). 

[5] M. Freitag, S. V. Kalinin, D. A. Bonnell, 
and A. T. Johnson, Phys. Rev. Lett. 89, 
216801 (2002). 

[6] C. Staii, A. T. Johnson, R. Shao, and 
D. A. Bonnell, NanoLett. 5, 893 (2005). 

[7] V. Meunier, S. V. Kalinin, J. Shin, 
A. P. Baddorf, and R. J. Harrison, Phys. 
Rev. Lett. 93, 246801 (2004); 

[8] Most commonly, 6 ~ 0.25 and I ~ 20 nm 
but much smaller values are feasible with 
ion-milled tips. 



[9] If the substrate has the dielectric constant 
e 1, its thickness is effectively reduced 
by e, e.g., for the geometrical thickness of 
200 nm and e = 3.9 of Si0 2 , D « 50 nm. 

[10] P. Pomorski, L. Pastewka, Ch. Roland, 
H. Guo, and J. Wang, Phys. Rev. B 69, 
115418 (2004). 

[11] M. Paillet, P. Poncharal, and A. Zahab, 
Phys. Rev. Lett. 94, 1868011 (2005), and 
references therein. 

[12] A. Svizhenko and M. P. Anantram, Phys. 
Rev. B 72, 85430 (2005). 

[13] J. D. Jackson, Am. J. Phys. 68, 789 (2000). 

[14] S. Belaidi, P. Girard, and G. Leveque, J. 
Appl. Phys. 81, 1023 (1997). 

[15] Single- Charge Tunneling, edited by 
H. Grabert and M. H. Devoret (Plenum, 
New York, 1992). 

[16] FastCap, Research Laboratory of 
Electronics, MIT. Available at 
http : //www . rle . mit . edu. 

[17] P. Jarillo-Herrero, S. Sapmaz, C. Dekker, 
L. P. Kouwenhoven, and H. S. J. van der 
Zant, Nature (London) 429, 389 (2004). 

[18] For a review, see S. M. Reimann and 
M. Manninen, Rev. Mod. Phys. 74, 1283 
(2002). 

[19] M. M. Fogler and E. Pivovarov, Phys. Rev. 
B 72, 195344 (2005). 

[20] K. A. Matveev, L. I. Glazman, and H. U. 
Baranger, Phys. Rev. B 54, 5637 (1996). 



9 



[21] W. G. van der Wiel, S. De Franceschi, J. 
M. Elzerman, T. Fujisawa, S. Tarucha, and 
L. P. Kouwenhoven, Rev. Mod. Phys. 75, 1 
(2003). 



10 



